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Summary. We present results from the first large-scale survey of neutron(n)- 
capture element abundances in planetary nebulae (PNe) . This survey was motivated 
by the fact that a PN may be enriched in n-capture elements if its progenitor star 
experienced s-process nucleosynthesis during the thermally-pulsing asymptotic gi- 
ant branch (AGB) phase. We have measured emission from Se and Kr in over 100 
PNe, and use literature data to expand our sample to 120 objects. [Kr III] 2.199 
and/or [Se IV] 2.287 /im were detected in 81 PNe, for a detection rate of nearly 70%. 
We derive Se and Kr abundances or upper limits using ionization correction factors 
derived from photoionization models. A significant range is found in the Se and Kr 
abundances, from near solar (no enrichment), to enriched by a factor of ten. Over- 
all, 41 of the 94 PNe with derived Se and/or Kr abundances or meaningful upper 
limits exhibit s-process enrichments. Our survey has increased the number of PNe 
with known ra-capture element abundances by an order of magnitude, enabling us 
to explore correlations between s-process enrichments and other nebular and central 
star properties. In particular, the Se and Kr enrichments display a positive correla- 
tion with nebular C/O ratios, as theoretically expected. Peimbert Type I PNe and 
bipolar PNe, whose progenitors are believed to be intermediate- mass stars (> 3- 

4 Mq), exhibit little or no s-process enrichment. Interestingly, PNe with H-deficient 
[WC] central stars do not exhibit systematically larger s-process enrichments than 
other PNe, despite the fact that their central stars are enriched in C and probably n- 
capture elements. Finally, the few PNe in our sample with known or probable binary 
central star systems exhibit little s-process enrichment, which may be explained if 
binary interactions truncated their AGB phases. We also briefly discuss a new ob- 
servational program to detect optical emission lines of n-capture elements, and new 
atomic data calculations that will greatly improve the accuracy of n-capture element 
abundance determinations in PNe. 
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1 Introduction 



Planetary nebulae (PNe) are the descendants of low- and intermediate-mass 
stars (M = 1-8 M Q ), the source of approximately half of the neutrons- 
capture element (atomic number Z > 30) nuclei in the Universe. These ele- 
ments are created by slow n-capture nucleosynthesis (the "s-process") during 
the thermally-pulsing asymptotic giant branch (AGB) phase, and are con- 
veyed to the stellar envelope during third dredge- up (TDU). Free neutrons 
are released during the intcrpulse phase by 13 C(a, n) 16 — or 22 Ne(a, n) 25 Mg 
in more massive AGB stars — and are captured by iron-peak "seed" nuclei, 
which undergo subsequent n-captures and /3-decays to transform into heav- 
ier elements [4, 10]. The enriched material is discharged into the ambient 
interstellar medium (ISM) via stellar winds and, ultimately, PN ejection. At 
solar metallicity, TDU occurs in stars with initial masses M > 1.5M ; in less 
massive stars, stellar winds reduce the mass of the stellar envelope below the 
critical mass for TDU (0.3-0.5 M ; [23]) before dredge-up can take place. 
Hence, stars less massive than 1.5 M Q (and their PNe) are not expected to 
exhibit s-process enrichments. 

Studying s-process enrichments in PNe provides unique information that 
cannot be obtained from spectroscopy of AGB stars. For example, the lightest 
n-capture elements (Z = 30-36) and noble gases generally cannot be detected 
in AGB stars, since strong lines of their neutrals and first ions reside in the 
UV where cool giant stars emit little flux. In addition, intermediate-mass stars 
(M > 3.5 Mq, hereafter IMS) are difficult to study during the AGB, due 
to heavy extinction from their dusty, optically thick circumstellar envelopes. 
However, these objects are readily observable as Type I PNe [16]. 



In this contribution, we present results from the first large-scale survey 
of n-capture elements in Galactic PNe. We have determined the elemental 
abundances (or upper limits) of the n-capture elements Se (Z = 34) and Kr 
(Z = 36) in 120 PNe, and examine correlations between s-process enrichments 
and other nebular and stellar properties for the first time. Se and Kr arc 
valuable tracers of s-process enrichments in PNe, since they are not expected 
to be depleted into dust; Kr is a noble gas, and Se has not been found to be 
depleted in the diffuse ISM [5] . A more detailed discussion of the results from 
our survey can be found in [21] and [22]. 

2 Observations and Abundance Determinations 

We have observed 103 Galactic PNe in the K band (2.14-2.30 fim) with 
the CoolSpec spectrometer [15] on the 2.7-m Harlan J. Smith telescope at 
McDonald Observatory. Including additional K band spectra from the lit- 
erature, our sample is comprised of 120 objects. We detected [Kr III] 2.199 
and/or [Se IV] 2.287 /xm in 81 of these objects, a remarkable detection rate 
considering the low cosmic abundances of Se and Kr (~ 2 x 10~ 9 relative to H 
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in the Solar System; [1]). These lines, first identified by [7], are resolved from 
other nebular features at our survey resolution of R — 500, with the exception 
of H 2 v = 3-2 lines in H 2 -emitting PNe (~30% of our targets). We removed 
the contaminating flux from these H 2 lines with the aid of high-resolution 
(R = 4400) observations and the measured fluxes of other observed H 2 lines 
[22]. 

Kr ++ and Se 3+ ionic abundances were derived for each object, using a 5- 
level and 2-level model atom, respectively. To determine elemental Se and Kr 
abundances, we utilized the ionization correction factors (ICFs) of [21], which 
were derived from Cloudy [8] and XSTAR [12] photoionization models. The 
derived Se and Kr abundances are accurate to within a factor of 2-3 for most 
objects in our sample, taking into account uncertainties in the line fluxes, T e , 
n e , and the fractional ionic abundances ++ /0, Ar ++ /Ar, and S ++ /S (taken 
from the literature) that arc incorporated into the ICFs. 

The choice of a reference element, for a metallicity- independent determi- 
nation of s-process enrichments, is particularly important. For most objects in 
our sample, we used O as a reference element since its abundance is generally 
more reliably determined than those of other elements. However, in Type I 
PNe, we found that [Ar/O] is larger by a factor of two relative to non-Type I 
PNe. Since Ar is not processed by PN progenitor stars, this implies that O 
was destroyed in these objects. O destruction can occur during hot bottom 
burning, if the temperature at the base of the convective envelope is high 
enough for the ON-cycle to be activated [13]. Therefore, we utilize Ar as a 
reference element for Type I objects, and use the notation [Se/(0, Ar)] and 
[Kr/ (O, Ar)] to underscore our choice of different reference elements for these 
two subclasses of PNe. 

Overall, we find a wide range of Se and Kr enrichment factors. The aver- 
age [Kr/(0, Ar)] = 0.98, with values for individual PNe ranging from —0.05 
to 1.89, while the mean [Se/(0, Ar)] = 0.31, with extrema at —0.56 and 
0.90 dex. In 18 PNe exhibiting both [Kr III] and [Se IV] emission, we find 
that [Kr/Se] = 0.5±0.2, in good agreement with theoretical predictions [3]. 
In 41 of the 94 PNe with derived abundances or meaningful upper limits, Se 
and/or Kr are enriched by more than a factor of two (the dispersion of light 
n-capture element abundances in unevolved solar-metallicity stars; [24]). We 
interpret the Se and Kr enrichments as evidence for in situ s-process nucle- 
osynthesis and TDU in PN progenitor stars. 



3 Correlations 

Our survey has increased the number of PNe with determined n-capture ele- 
ment abundances by nearly a factor of ten. This enables us to search for corre- 
lations between s-process enrichments and other nebular and stellar properties 
for the first time. In Table 1, we show the mean values of [Se/(0, Ar)] and 
[Kr/(0, Ar)] for different subclasses of PNe. 
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Table 1. Mean Se and Kr Abundances in Different Subclasses of PNe a 





Mean 




Number of 


Mean 




Number of 


Property 


[Se/(0, Ar)] < o > b 


Se Detections [Kr/(0, Ar)] < o > b 


Kr Detections 


Progenitor Mass 


Type I 


-0.03 


0.27 


12 


0.09 


0.14 


3 


Non-Type I 


0.36 


0.26 


55 


1.02 


0.27 


30 


Morphology 


Bipolar 


0.27 


0.38 


14 


0.68 


0.25 


8 


Elliptical 


0.28 


0.22 


28 


1.09 


0.38 


15 


Central Star Type 


[WC] 


0.39 


0.28 


16 


0.90 


0.34 


10 


Non-[WC] 


0.29 


0.28 


39 


1.05 


0.32 


20 


Binary 


-0.07 


0.12 


5 


0.82 


0.08 


2 


Full Sample 


0.31 


0.27 


67 


0.98 


0.31 


33 



(a) Only PNe exhibiting Se and/or Kr emission and with determined O and Ar 
abundances are considered. 

(b) The < a > are mean absolute deviations in the Se and Kr abundances. 



Theoretically, it is expected that s-process enrichments will be correlated 
with the C/O ratio, since 12 C is transported to the stellar envelope along 
with n-capture elements during TDU [4] . Strong empirical evidence has been 
found for such a correlation in AGB [20] and post-AGB stars [25]. Indeed, 
[Se/(0, Ar)] and [Kr/(0, Ar)] are positively correlated with the gaseous C/O 
ratio (Figure 1), with correlation coefficients r = 0.45 and 0.64, respectively. 
Interestingly, Se and Kr enrichments do not increase with C/O as rapidly as 
Sr, Y, and Zr (observed in AGB stars). This is likely an effect of the smaller 
s-process yields of Se and Kr relative to these three elements [3]. 

We also find that Type I PNe, which are descendants of IMS, exhibit little 
s-process enrichment compared to non-Type I PNe (Table 1, Figure 2). Bipo- 
lar PNe also show smaller s-process enrichments than elliptical PNe, although 
the discrepancy is not as pronounced in this case. This has been verified with 
Kolmogorov-Smirnov (KS) tests, which show that the probabilities that the 
Se and Kr enrichments of Type I and non-Type I PNe are drawn from the 
same cumulative distribution functions are p^ s = 0.02 and p^s = 0.01, respec- 
tively; for bipolar and elliptical PNe, pks(Se) = 0.42 and pks(Kr) = 0.21. The 
low enrichment factors for Type I and bipolar PNe are likely due to the small 
masses of their progenitor stars' intershell layers (relative to those of lower 
mass AGB stars), as well as the severe dilution that the enriched material 
experiences when it is dredged-up into the massive envelopes of these stars 
[14]. This result is in agreement with the marginal Zr enrichments found by 
[9] in a sample of intermediate-mass AGB stars. 

[WC] PN central stars are H-deficient objects whose surface compositions 
are enriched in C and (likely) n-capture elements [27]. It may therefore be 




Fig. 1. [Se/(0, Ar)] (left) and [Kr/(0, Ar)] (right) plotted against the gaseous 
C/O ratio of PNe in our sample. The best linear fit to each correlation is shown as a 
solid line (with correlation coefficient r and significance p r =o indicated) . The dashed 
lines are fits to [<Sr, Y, Zr>/Fe] as a function of C/O in AGB and post-AGB stars 
[19, 20, 18, 26]. 



expected that their surrounding nebulae are also enriched in these nuclei. 
However, we find that there is no systematic difference between the Se and 
Kr enrichments of [WC] PNe and nebulae with H-rich central stars (Table 1): 
KS tests show that there is an 85% and 99% probability that the Se and 
Kr abundances (respectively) in these two PN subclasses are drawn from the 
same cumulative distribution functions. It has previously been shown that 
the compositions of [WC] PNe are similar to PNe with H-rich nuclei for other 
elements [17], even C [6]. 

Relatively few PNe in our sample exhibit direct evidence for binary central 
star systems. Izzard ([11], these proceedings) showed that binary interactions 
that take place during the thermally-pulsing AGB can truncate this phase of 
evolution by enhancing the mass-loss of the AGB star. If this occurs, these 
systems will be less enriched in TDU products than otherwise similar sin- 
gle stars. Interestingly, the 14 objects in our sample that display evidence of 
binary central stars exhibit little to no s-process enrichment (Table 1). How- 
ever, this is not necessarily due to the binary companions, since the level of 
enrichment is also strongly dependent on the PN progenitor's initial mass. 



4 Conclusions and Future Work 



We have highlighted the major results from our large-scale survey of near- 
infrared [Kr III] and [Se IV] emission lines in Galactic PNe. We detected Se 
and/or Kr emission in 81 out of 120 objects in our sample, of which 41 are 
significantly enriched due to in situ s-process nucleosynthesis and TDU in 
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Full Sample 



□ with Type I/Bipolar PNe 
H Without Type 1/ 
Bipolar PNe 



Full Sample 





Fig. 2. Histograms of Se and Kr abundances, separated into 0.1 dex bins, are 
shown for the full sample (top panels), Type I PNe (middle panels), and bipolar 
PNe (bottom panels). 



their progenitor stars. Se and Kr abundances are correlated with the gaseous 
C/O ratio, as predicted by theoretical nucleosynthetic models. Type I and 
bipolar PNe, whose progenitors are IMS, exhibit little if any s-process enrich- 
ment compared to objects with less massive progenitors. On the other hand, 
we find that [WC] PNe are not systematically more s-process enriched than 
objects with H-rich nuclei, which stands in stark contrast to the strong C 
and probable s-process enrichments of [WC] central stars. PNe with binary 
central star systems show little evidence of s-process enrichments, as may be 
expected if binary interactions truncated their thermally-pulsing AGB phase. 

The derived Se and Kr abundances are uncertain by a factor of 2-3 for most 
objects in our sample. We have shown that these uncertainties arise primarily 
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from the ICFs [21], for two reasons. First, we have detected only one ion 
each of Se and Kr. Therefore, the ICFs can be large and the uncertainties 
significant. Secondly, the atomic data controlling the ionization balance of Se 
and Kr — photoionization (PI) cross-sections and rate coefficients for various 
recombination processes — are poorly if at all known (indeed, this is true 
of most n-capture element ions). To derive the Se and Kr ICFs, we used 
approximations for these atomic data in Cloudy and XSTAR [21]. 

We have recently instigated a new project to improve the accuracy of 
n-capture element abundance determinations in PNe. We are observing s- 
process enriched PNe drawn from our near-infrared sample in the optical, 
in order to search for lines from additional ions of Kr and Se (as well as 
transitions of other n-capture elements). In the optical spectra of five PNe, 
we have detected [Kr IV] in four and [Kr V] in one thus far. The detec- 
tion of additional ions reduces the magnitude and importance of uncertainties 
in the ICFs. In addition, we are computing new atomic data for the first 
six ions of Se and Kr with the atomic structure code AUTOSTRUCTURE 
[2], including PI cross-sections and rate coefficients for radiative and dielcc- 
tronic recombination. These calculations are complemented by experimental 
absolute PI cross-section measurements near the ionization threshold of these 
ions, performed at the Advanced Light Source synchrotron radiation facility 
in Berkeley, CA. The new atomic data determinations will enable us to de- 
rive more accurate and robust ICFs with photoionization models. Combined 
with the new observational data, this will significantly reduce uncertainties in 
Se and Kr abundance determinations in PNe and allow for a more rigorous 
investigation of s-process enrichments. 

This work has been supported by NSF grants AST 97-31156 and AST 04- 
06809. 
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